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Abstract 

We study v e and v e signals from a supernova with strong magnetic field detected by a 100 kton liquid Ar detector. The 
change of neutrino flavors by resonant spin-flavor conversions, matter effects, and neutrino self-interactions are taken 
into account. Different neutrino signals, characterized by the total event numbers of v e and v e and neutronization burst 
event, are expected with different neutrino oscillation parameters and neutrino magnetic moment. Observations of 
supernova neutrino signals by a 100 kton liquid Ar detector would constrain oscillation parameters as well as neutrino 
magnetic moment in either normal and inverted mass hierarchies. 
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1. Introduction 

Identification of very large neutrino magnetic moment (/i v » 10~ 18 jUb, where fig is Bohr magnetons) is one of 
important traces of particle physics beyond the standard model. Neutrino experiments have constrained the upper 
limit of the neutrino magnetic moment as fi Ve < 5.8 x 10~ U jUb I ID- If neutrinos are Majorana particles and they have 
magnetic moment, they have only transition magnetic moment. The magnetic moment enables a spin precession 
between a left-handed neutrino Vl and a right-handed antineutrino vr with different flavors in a strong magnetic field 
(e.g., J2l)- This is the resonant spin-flavor (RSF) conversion. 

Supernova (SN) neutrino is a good indicator observing the effects of neutrino magnetic moment. The magnetic 
field of Fe core of a presupernova is expected to attain 10 1() gauss yfl. Proto-neutron stars should have stronger mag- 
netic field since observations of pulsars have magnetic field of ~ 10 12 gauss. If neutrinos are Majorana particles and 
have large transition magnetic moment (ju v ~ I0^ n fig), the RSF conversion is expected to occur in innermost region 
of SN ejecta (e.g., no, |g, |7[]). The final neutrino spectra also have influences of Mikheyev-Smirnov-Wolfenstein 
(MSW) effects (e.g., lU) and neutrino self-interactions (e.g., |@]). The influence of the RSF conversion, which will be 
observed in water-Cerenkov detector, on SN v e signal has been predicted flBl l6l [Toll . The observations of v e signals 
are expected to reveal the effects of the RSF conversions in inverted mass hierarchy. 

A liquid Ar detector distinguishes v e and v e events by charged-current (CC) reactions, events by neutral-current 
(NC) reactions and electron scattering of all flavors of neutrinos and antineutrinos. Especially, many events of v e 
will be observed owing to large cross section of 40 Ar(v e , e~) reaction. Previous studies have estimated the SN neutrino 
events by liquid Ar detector taking into account the MSW effects fill [Till and neutrino self-interactions lfl3ll . ICARUS 
experiment (e.g., ifUl l is proposed as the usage of a 3 kton liquid Ar detector. The experiment with the 600 ton detector 
T600 started at the underground laboratory in Gran Sasso 11511 . In Long-Baseline Neutrino Experiment Project it is 
proposed to construct three 1 7-kton liquid Ar detectors at DUSEL lfl6Tl . Development of a very large liquid Ar detector 
with 100 kton is expected. 

We investigate possible SN neutrino signals observed by a 100 kton liquid Ar detector taking into account the RSF 
conversions, MSW effects, and neutrino self-interactions. We evaluate the SN neutrino events for v e and v e by the 
CC reactions and all flavors of neutrinos and antineutrinos by the NC reactions. In this Letter, we report for the first 
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Figure 1 : Schematic diagram of effective masses of neutrinos and antineutrinos in (a) normal mass hierarchy and (b) inverted mass hierarchy. Solid 
circles indicate three RSF resonances (RSF-H, RSF-X, and RSF-L) and dashed circles indicate two MSW resonances (MSW-H and MSW-L). The 
changes of effective mass in accordance with a solid curve and dashed curve indicate adiabatic and non-adiabatic conversions, respectively. 



time the possibility that the effects of the RSF conversions and large neutrino magnetic moment are extracted from 
SN signals in both normal and inverted mass hierarchies. 



2. Model 

2.1. SN neutrino model 

We used a SN model adopted in 0ll[l2l]- The SN progenitor is evolved from a 15 M Q star OH. We used the SN 
density profiles at 0, 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10 s after the core bounce IU7I1 to see the shock effect of the RSF 
conversions and the MSW effects. When the RSF conversions are taken into account, the SN magnetic field transverse 
to the neutrino propagation at the radius r was assumed to be Bj_(r) = 10 n (10 8 cm/r) 3 gauss llllfl]. The electron 
fraction in the innermost hot-bubble and wind region where the mass coordinate is smaller than 1 .43 M was assumed 
to be Y e = 0.49. The electron fraction of the outer region was assumed to be the same as that of the progenitor model. 
We assumed that the SN magnetic field and the electron fraction do not change after the shock arrival. We should note 
that the magnetic field and the electron fraction after the shock propagation are still uncertain. They would change 
complicatedly. The distance from the SN was set to be 10 kpc. 

We assumed the time evolution of the neutrinos emitted from the neutrino sphere as follows. The neutrino flux 
decreases exponentially with the time scale of t v — 3 s. The total neutrino energy is set to be 3 x 10 53 ergs, correspond- 
ing to the binding energy of a 1.4 M neutron star (e.g., ifioll ). The neutrino energy is equipartitioned to each flavor 
of neutrinos. We considered two spectral types of the neutrino energy spectra, which obey Fermi distributions. The 
neutrino temperatures in Spectral Type 1 were set to be (T Ve , T Ve , T v ) - (3.2 MeV, 5 MeV, 6 MeV), where T Ve , T Ve , 
and T Y mean the temperatures of v e , v e , and v MiT and v^ T , respectively |T(1 19]. The neutrino chemical potential of 



each flavor was set to be zero. In Spectral Type 2, the temperatures and chemical potentials of the neutrinos were set 
as (T Ve , T Ve , r Vr ) = (3.5 MeV, 4 MeV, 7 MeV) and (/^//y^/V) = (7.4 MeV, 10 MeV, MeV). This set corresponds 
to the neutrino spectra of the Livermore model 1171 12011 . 



In the case of inverted mass hierarchy, neutrino self-interactions change neutrino-flavors in the innermost region 
inside RSF-H resonance. We considered the flavor conversion by neutrino interactions in inverted mass hierarchy. 



The energy splitting in neutrino spectra was taken into account in accordance with [21]. Electron neutrinos above the 
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Table 1: The final neutrino flux of v e [0 Vf (fin)] and v e [<f>y e (fin)] with the relation of the original neutrino flux v a [<p Va ]. The neutrino flux <py x 
means the original flux of Vu, v T , v„, v r . The first column is model of neutrino oscillation parameters: the adiabaticity of the RSF-H resonance (R 
or M), the adiabaticity of the MSW-H resonance (A or N), and the mass hierarchy (N or I). For example, RAN means adiabatic RSF-H resonance, 
adiabatic MSW-H resonance, and the normal mass hierarchy. In inverted mass hierarchy, the v e flux above the critical energy Ec is written. See 
text for details. 
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critical energy Ec are converted from v^ T . In the case of anti-neutrinos, v e is converted from v MiT approximately in all 
energy range. The critical energy Ec is 8. 16 MeV and 7.55 MeV for Spectral Types 1 and 2, respectively. 



We set the neutrino oscillation parameters in accordance with the observational constraints (Hill IH. The 
squared mass differences Am?. = mj - m 2 . were set to be Am 2 j = 7.6 x 10~ 5 eV 2 and |Am 2 j| = 2.5 x 1CT 3 eV 2 . 

The mixing angles 9\o and 023 were set to be sin 2 20i2 = 0.87 and sin 2 2023 = 1. The mixing angle 0^ was set 
to be sin 2 20i3 = 0.04 as a large value corresponding to an adiabatic MSW-H resonance and sin 2 20i3 = 10~ 6 as a 
small value corresponding to a non-adiabatic MSW-H resonance. Mass hierarchy was also considered as a parameter. 
When the RSF conversion was taken into account, the strength of the neutrino transition magnetic moment was set to 
be 10" 12 /^b- We numerically solved the neutrino flavor change in accordance with ifioll . 

2.2. Neutrino detection 

We considered the neutrino detection by a liquid Ar detector. The fidual volume of the liquid Ar detector was 
assumed to be 100 kton. This is about 30 times larger than the volume of ICARUS liquid Ar detector. The detection 
efficiency of electrons and positrons is assumed to be unity when the energy of electrons and positrons is larger than 



E t h = 5 MeV and otherwise zero. The energy resolution was taken as A(E e ) = 0.11 V£' e (MeV) + 0.025£' c MeV 
In this study, we evaluated the total event numbers of the CC reactions °Ar(v e , e~), 40 Ar(v f , e + ), and the NC reactions 
40 Ar(v, v'). The rates of the neutrino events and electron and positron energy spectra of them were evaluated using the 
same procedure described in [10]. The cross sections of the neutrino reactions with 40 Ar were taken from | IT ] . 



3. Results 

Flavor changes by the RSF conversions and the MSW effects are very complicated. The neutrino flavors are 
converted in accordance with schematic diagram of effective masses of neutrinos and antineutrinos in Fig. 1 (see also 
Si) . This figure also indicates the locations of three RSF resonances, i.e., RSF-H, RSF-X, and RSF-L (e.g., fiSH) 
and two MSW resonances, i.e., MSW-H and MSW-L (e.g., 181]). The adiabaticities of these resonances depend on the 
electron-number density gradient, neutrino energy, the mixing angles, the mass difference, and the stellar magnetic 
field (e.g., J10S]). 

We evaluated the final flux of v e and v e as a function of the initial neutrino flux. The relation of the final flux of 
v e and v e to the initial neutrino flux is listed in Table 1. Here we assumed that \U e \\ 2 — cos 2 0i2, |£/ e 2l 2 = sin 2 0i2, and 
It^rfl 2 = 0, respectively, for simplicity. We also assumed that the resonances of RSF-X and RSF-L are non-adiabatic 
and the MSW-L resonance is adiabatic. We note in inverted mass hierarchy that the neutrino flux above the critical 
energy Ec is converted as (0 V ,, <f> v , Vr ) — > ((0 V +0vJ/2, (0 V , +0 V(i )/2, (0 V , + Vr )/2) and antineutrino flux is converted 
as (0p,,0v , Vr ) — > ((0i>„ +0v r )/2, (0v, +0i>„)/2, (0i>, + 0v T )/2)by neutrino self-interactions before the RSF conversions. 
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Table 2: The event numbers of CC 40 Ar(v e , e ), CC 40 Ar(v c , e + ), and NC 40 Ar(v, v') reactions with the relation to models of neutrino oscillation 
parameters. 



Reaction RAN RAI RNN RNI MAN MAI MNA MNI 

Spectral Type 1: (J Ve , T Ve , T V/it ) =(3.2 MeV, 5 MeV, 6 MeV), (ji Ve , fi Ve , f^J =(0 MeV, MeV, MeV) 
CC 40 Ar(v e , O 16384 '' 11970 12376 11981 16391 12585 12312 12541 
CC 40 Ar(v e ,e + ) 762 763 766 819 580 672 578 823 

NC 40 Ar(v,vQ 13668 13734 13663 13733 13805 13838 13805 13805 

Spectral Type 2: (J v „ T- v „ T V J =(3.5 MeV, 4 MeV, 7 MeV), Qj v „ =(7.4 MeV, 10 MeV, MeV) 

CC 4U Ar(v e ,e-) 19157 15534 16656 13534 21908 16950 16579 16981 
CC 40 Ar(v e , e + ) 1105 1119 1114 1289 603 866 596 1294 
NC 40 Ar(v,v') 19063 19255 19047 19253 19463 19464 19464 19463 



In Table 1, there are eight models of neutrino oscillation parameters: the adiabaticities of the RSF-H and MSW- 
H resonances and the mass hierarchies. We denote these models as three characters. The first character means the 
adiabaticity of the RSF-H resonance; "R" means that the RSF-H is adiabatic, i.e., the RSF conversion is effective, 
and "M" means that the RSF-H resonance is non-adiabatic thus the RSF conversion does not work. The second 
character means the adiabaticity of the MSW-H resonance; "A" means adiabatic and "N" means non-adiabatic. The 
third character means the mass hierarchy; "N" and "I" correspond to "normal" and "inverted" mass hierarchies, 
respectively. 

Generally, the temperatures of SN neutrinos indicate the relation of T Ve < Ty e < T v r . In this case the event 
number of v e , N Ye , or v e , Ny r , becomes large when v e or v e is converted from energetic v^ T . From Table 1, we expect 
the event number of v e increases with [pattern a] JV Vf (RAI), JV Vf (RNI) < [pattern b] N Ve (MAI), jV Ve (MNI) < [pattern c] 
Af Ve (RNN), jV Vf (MNN) < [pattern d] N Ve (RAN) < [pattern e] N Ve (MAN) and that the event number of v e increases with 
[pattern a] jvJ(MAN), JVy c (MNN) < [pattern b] M> C (MAI) < [pattern c] M> C (RAI), iVy f (RNI), M> e (MNI) < [pattern d] 
Ny e (RAN), A^ Ve (RNN). 

Considering the above relations we numerically evaluated the total event numbers of the CC reactions 40 Ar(v c , e~) 
and 40 Ar(v e , e + ) and the NC reactions 40 Ar(v, V) with the two SN neutrino spectral types. Table 2 shows the total 
event numbers of these three reactions. The event numbers of 40 Ar(v e , e~) and the NC reactions are of order 10 4 . The 
event number of 40 Ar(v e , e + ) is in the range of hundreds to 1300. The event numbers of the CC reactions depend on 
oscillation parameters and neutrino magnetic moment. However, the NC event scarcely depends on. When a model 
of oscillation parameters is given, the event numbers in Spectral Type 2 are larger than the corresponding numbers 
in Spectral Type 1. This relation is due to the fact that the average energies of v e , v e , and v^ T in Spectral Type 2 are 
larger than the corresponding average energies in Spectral Type 1 . 

Next, we indicate the dependence of the event number ratio of 40 Ar(v c , e~) reaction to the NC reactions (rcc(? e )/nc = 
NvJNnc) and that of 40 Ar(v e , e + ) reaction to the NC reactions (/ccoy/NC = NyJN^c) on the oscillation parameters. 
Figure 2(a) shows the relation of the two event number ratios in Spectral Type 1 . The event number ratios on this 
figure do not overlap with two exceptions. The ratios of models RNN and MNI overlap with error bars. The ratios of 
models RAI and RNI also do. 

The event number ratios in Fig. 2(a) do not obey the inequalities indicated above. This is due to the shock effect 
on the RSF conversions. When the shock wave passes through an adiabatic resonance, the adiabaticities changes to 
non-adiabatic. Since the RSF-H resonance is located in deep region of the SN ejecta, the adiabaticity of the resonance 
changes when neutrino irradiation is still strong. As a result, the shock effect at the RSF-H resonance is stronger than 
that at the MSW-H resonance in the total neutrino event. 

We also show the event ratios as open symbols in models RAN, RAI, RNA, and RNI without adiabaticity change 
by the shock passage of the RSF-H resonance. The ratios with closed symbols in the four models are different from 
the corresponding ratios with relations of the event numbers estimated in Table 1. When the shock effect is not 
taken into account, the numerical result is consistent with the theoretical estimation. For example, the event ratio 
for v e increases with [pattern a] JV Ve (RAI), JV Ve (RNI) < [pattern b] JV Ve (MAI), AT Ve (MNI) and [pattern c] JV Vf (RNN), 
Af V( ,(MNN) < [pattern d] N Vf (RAN) < [pattern e] jV Vf (MAN). 
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Figure 2: The relation between the event number ratio of 40 Ar(v,,, e~) to the NC reactions of 40 Ar to the NC reactions of 40 Ar rcc(v e )/NC and 'he 
ratio of 40 Ar(V(,, e + ) to the NC reactions of 40 Ar f"cc(i> e )/NC- Panels (a) and (b) are the number ratios in Spectral Types 1 and 2, respectively. Squares, 
circles, and triangles indicate patterns A, B, and C, respectively, of the observation of neutronization burst (see text for more details). Open symbols 
indicate the corresponding ratios when the change of the adiabaticities in the RSF-H and MSW-H resonances by the shock propagation in the SN 
is not taken into account. 



When the RSF conversions are effective, the event ratios with shock effect are different from those without shock 
effect. On the other hand, the ratios with ineffective RSF conversion cases hardly depend on the shock effect. This is 
because the time when the shock arrives at the MSW-H resonance is later than the shock arrival time to the RSF-H 
resonance and most of neutrinos have passed until the shock arrival at the MSW-H resonance. If we obtain large 
difference between the event ratios for v e and v e with the shock effect and those without the shock effect, we can 
discuss the possibility of the RSF conversions. 

Neutronization burst is another special neutrino signal appearing at the core bounce. The spectra of the neutron- 
ization burst are quite different from the total neutrino signal and the main component of the neutronization burst is 
v e . The neutrinos in the neutronization burst are scarcely affected by neutrino self-interactions. Most of the released 
neutrinos are v e in neutronization burst, so that the interaction of v e and other flavor of neutrinos and antineutrinos 
is quite small. The neutrino signal depends on neutrino oscillation parameters. The effects of RSF conversions and 
the MSW effects on v e have been investigated in Q4J, . We investigate both v e and v e signals in the neutronization 
burst. Figure 2 shows three patterns of the detection of the neutronization burst. Table 3 indicates the final flux of 
v e and v e . Pattern A corresponds to no signals of neutronization burst. Pattern B corresponds to the v e detection 
and Pattern C indicates v e detection in neutronization burst. The signal of the neutronization burst distinguishes the 
neutrino events between model RAI and model RNI. Thus, the v e and v e signals from the neutronization burst also 
constrain the neutrino oscillation parameters and neutrino magnetic moment. 

Figure 2(b) shows the relations of the two event number ratios in Spectral Type 2. We see that the relation of the 
event number ratios is similar to Fig. 2(a); the dependence of the event number ratios on oscillation parameters is 
larger. This is due to larger difference between the v UiT temperature and those of v e and v e . 

If we can expect the initial neutrino spectra precisely, we would constrain the neutrino oscillation parameters and 
neutrino magnetic moment from SN neutrino signals detected by a liquid Ar detector. The total events of the CC 
reactions of v e and v e and the NC reactions will indicate different signals by different oscillation parameters for most 
cases. The shock effect in the total CC events of v e and v e will clarify whether the RSF conversions are effective or not. 
This suggests the upper limit of the product of neutrino magnetic moment and the SN magnetic field. The detection of 
neutronization burst signal by v e or v e detector or no detection gives additional constraint to the oscillation parameters. 

Even if the initial neutrino spectra are expected less precisely, we would constrain some oscillation parameters and 
neutrino magnetic moment. Large and small v e events with large v e events and no detection of neutronization burst 
should be the evidence for models RAN and MAN, respectively. Small v e event, large v e event, and no detection of 
neutronization burst will indicate the evidence for model RNI. As first discussed in [6], the detection of neutronization 
burst in v e signal is a strong evidence for model RAI. Small v e event will confirm this evidence. Small v e event with 
the detection of the neutronization burst in v e signal will exclude the possibility of large neutrino magnetic moment 
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Table 3: The final neutrino flux of v e [0,,,(fin)] and v e [0,.,(fin)] with the relation of the original neutrino flux v a [<p Va ] in the neutronization burst. 
The first column is model of neutrino oscillation parameters as described in Table 1 . 
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enough to bring about the RSF conversions. 

We have shown that the shock effect appears in the event number ratios of v e and v e signals when the RSF 
conversions are effective. However, if a SN explodes jet-like, different time dependence will be seen by the adiabaticity 
change of the MSW-H resonance B25I1 . Magneto-hydrodynamic simulations of SNe indicated that magnetically-driven 
jet explosion is induced by strong toroidal magnetic field [26]. Under such explosions the points and adiabaticities 
of the RSF and MSW resonances would be quite different from those in spherical explosions. Thus, the investigation 
of the neutrino signal from magnetically-driven SNe taking into account the RSF conversion is fascinated. Analyses 
with optical observations of SNe 12711 help distinguishing the effects by the RSF conversions and jet-like explosion of 
theSN. 



4. Conclusions 

We have investigated traces of large neutrino magnetic moment, the RSF conversions, in SN neutrino signals 
observed by a 100 kton liquid Ar detector. We predict that different signals will be observed in the event number 
ratios of the CC v e and v e reactions to the NC reactions with different oscillation parameters in both normal and 
inverted mass hierarchies. Large shock effect will be found in the total event number ratios when the RSF conversions 
are effective. Combined with the signals of the neutronization burst, the v e and v e signals would constrain neutrino 
oscillation parameters and the neutrino magnetic moment. The observations of both v e and v e will open the possibility 
for constraining neutrino oscillation parameters as well as neutrino magnetic moment. 

T.Y. thanks Mary Bishai and Kate Scholberg for kindly providing useful information of liquid Ar detectors in 
LBNE Project. Numerical computations were in part carried out on the general-purpose PC farm at Center for Com- 
putational Astrophysics, CfCA, of National Astronomical Observatory of Japan. This work has been supported in part 
by the Ministry of Education, Culture, Sports, Science and Technology, Scientific Research (C)(20540284, 23540287). 



References 

[1] A. G. Beda, et al., Phys. At. Nucl. 70 (2007) 1873. 

[2] C. S. Lim, W. J. Marciano, Phys. Rev. D 37 (1988) 1368. 

[3] A. Heger, S. E. Woosley, H. C. Spruit, Astrophys. J. 626 (2005) 350. 

[4] S. Ando, K. Sato, Phys. Rev. D 67 (2003) 023004. 

[5] S. Ando, K. Sato, Phys. Rev. D 68 (2003) 023003. 

[6] S. Ando, K. Sato, J. Cosmo. Astropart. Phys. 10 (2003) 001. 

[7] E. K. Akhmedov, T. Fukuyama, J. Cosmo. Astropart. Phys. 12 (2003) 010. 

[8] A. S. Dighe, A. Y. Smirnov, Phys. Rev. D 62 (2000) 033007. 

[9] Y.-Z. Qian, G. M. Fuller, Phys. Rev. D 51 (1995) 1479. 

[10] T. Yoshida, A. Takamura, K. Kimura, H. Yokomakura, S. Kawagoe, T. Kajino, Phy. Rev. D 80 (2009) 125032. 

[11] I. Gil-Botella, A. Rubbia, J. Cosmology Astropart. Phys. 10 (2003) 009. 

[12] I. Gil-Botella, A. Rubbia, J. Cosmology Astropart. Phys. 08 (2004) 001. 

[13] S. Choubey, B. Dasgupta, A. Dighe, A. Mirizzi, arXiv: 1008.0308 (2010). 



6 



[14] The ICARUS-Milano Collaboration, F. Arneodo, et al., Eur. Phys. J. C 33 (2004) 233. 

[15] The ICARUS-Milano Collaboration, F. Arneodo, et al., Phys. Rev. D 74 (2006) 1 12001. 

[16] M. Bass et al., The Interim LBNE Physics Report (2010). 

[17] S. Kawagoe, T. Yoshida, T. Kajino, H. Suzuki, K. Sumiyoshi, S. Yamada, Phys. Rev. D 81 (2010) 123014. 

[18] S. E. Woosley, T. A. Weaver, Astrophys. J. Suppl. 101 (1995) 181. 

[19] T. Yoshida, T. Kajino, H. Yokomakura, K. Kimura, A. Takamura, D. H. Hartmann, Phys. Rev. Lett. 96 (2006) 091 101. 

[20] T. Totani, K. Sato, H. E. Dalged, J. R. Wilson, Astrophys. J. 496 (1998) 216. 

[21] G. L. Fogli, E. Lisi, A. Marrone, A. Mirizzi, J. Cosmology Astropart. Phys. 12 (2007) 010. 

[22] S. N. Ahmed et al. (SNO Collaboration), Phys. Rev. Lett. 92 (2004) 181301. 

[23] J. Hosaka et al. (Super-Kamiokande Collaboration), Phys. Rev. D 74 (2006) 032002. 

[24] S. Abe et al. (KamLAND Collaboration), Phys. Rev. Lett. 100 (2008) 221803. 

[25] S. Kawagoe, T. Takiwaki, K. Kotake, J. Cosmology Astropart. Phys. 09 (2009) 033. 

[26] T. Takiwaki, K. Kotake, K. Sato, Astrophys. J. 691 (2009) 1360. 

[27] K. Maeda, et al., Science 319 (2008) 1220. 



7 



